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The production of jet pairs at large rapidity difference at hadron colliders is potentially sensitive to BFKL physics. We 
present the results of a BFKL Monte Carlo calculation of dijets at the Tevatron. The Monte Carlo incorporates kinematic 
effects that are absent in analytic BFKL calculations; these effects significantly modify the behavior of dijet cross sections. 



1. MONTE CARLO APPROACH TO BFKL 

Fixed-order QCD perturbation theory fails in some 
asymptotic regimes where large logarithms multiply 
the coupling constant. In those regimes resummation 
of the perturbation series to all orders is necessary to 
describe many high-energy processes. The Balitsky- 
Fadin-Kuraev-Lipatov (BFKL) equation [ [lj] performs 
such a resummation for virtual and real soft gluon 
emissions in such processes as dijet production at large 
rapidity difference in hadron-hadron collisions. BFKL 
resummation gives [ ^[ a subprocess cross section that 
increases with rapidity difference as a ~ exp(AA), 
where A is the rapidity difference of the two jets with 
comparable transverse momenta pri and PT2- 

Experimental studies of these processes have re- 
cently begun at the Tevatron pp and HERA ep col- 
liders. Tests so far have been inconclusive; the data 
tend to lie between fixed-order QCD and analytic 
BFKL predictions. However the applicability of an- 
alytic BFKL solutions is limited by the fact that they 
implicitly contain integrations over arbitrary numbers 
of emitted giuons with arbitrarily large transverse mo- 
mentum: there are no kinematic constraints included. 
Furthermore, the implicit sum over emitted giuons 
leaves only leading-order kinematics, including only 
the momenta of the 'external' particles. The absence 
of kinematic constraints and energy-momentum con- 
servation cannot, of course, be reproduced in exper- 
iments. While the effects of such constraints are in 
principle sub-leading, in fact they can be substantial 
and should be included in predictions to be compared 
with experimental results. 

The solution is to unfold the implicit sum over glu- 
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ons and to implement the result in a Monte Carlo 
event generator [ ||, [|. This is achieved as follows. 
The BFKL equation contains separate integrals over 
real and virtual emitted giuons. We can reorganize 
the equation by combining the 'unresolved' real emis- 
sions — those with transverse momenta below some 
minimum value (chosen to be small compared to the 
momentum threshold for measured jets) — with the 
virtual emissions. Schematically, we have 



virtual 



real 



virtual+real ,unres . 



(1) 



real,res. 



We perform the integration over virtual and unresolved 
real emissions analytically. The integral containing the 
resolvable real emissions is left explicit. 

We then solve by iteration, and we obtain a differ- 
ential cross section that contains a sum over emitted 
giuons along with the appropriate phase space factors. 
In addition, we obtain an overall form factor due to 
virtual and unresolved emissions. The subprocess cross 
section is 



da = d&o x \ f n 



(2) 



where /„ is the iterated solution for n real giuons emit- 
ted and contains the overall form factor. It is then 
straightforward to implement the result in a Monte 
Carlo event generator. Because emitted real (resolved) 
giuons appear explicitly, conservation of momentum 
and energy, as well as evaluation of parton distribu- 
tions, is based on exact kinematics for each event. In 
addition, we include the running of the strong coupling 
constant. See [ || for further details. 

2. DIJET PRODUCTION AT HADRON 
COLLIDERS 

At hadron colliders, the BFKL increase in the dijet 
subprocess cross section with rapidity difference A is 
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Figure 1. The azimuthal angle decorrelation in dijet production at the Tevatron as a function of dijet rapidity 
difference A, for jet transverse momentum px > 20 GeV. The analytic BFKL solution is shown as a solid curve 
and a preliminary D0 measurement [ |(| is shown as diamonds. Error bars represent statistical and uncorrelated 
systematic errors; correlated jet energy scale systematics are shown as an error band. 



unfortunately washed out by the falling parton distri- 
bution functions (pdfs). As a result, the BFKL pre- 
diction for the total cross section is simply a less steep 
falloff than obtained in fixed-order QCD, and tests of 
this prediction are sensitive to pdf uncertainties. A 
more robust pediction is obtained by noting that the 
emitted gluons give rise to a decorrelation in azimuth 
between the two leading jets.[ ||, || This decorrelation 
becomes stronger as A increases and more gluons are 
emitted. In lowest order in QCD, in contrast, the jets 
are back-to-back in azimuth and the (subprocess) cross 
section is constant, independent of A. 

This azimuthal decorrelation is illustrated in Fig- 
ure |] for dijet production at the Tevatron pp collider [ 
D, with center of mass energy 1.8 TeV and jet trans- 
verse momentum px > 20 GeV. The azimuthal an- 
gle difference A(j> is defined such that cosA^ = 1 for 
back-to-back jets. The solid line shows the analytic 
BFK1 prediction. The BFKL Monte Carlo prediction 
is shown as crosses. We see that the kinematic con- 



straints result in a weaker decorrelation due to sup- 
pression of emitted gluons, and we obtain improved 
agreement with preliminary measurements by the D0 
collaboration [^], shown as diamonds in the figure. 

In addition to studying the azimuthal decorrelation, 
one can look for the BFKL rise in dijet cross section 
with rapidity difference by considering ratios of cross 
sections at different center of mass energies at fixed 
A. The idea is to cancel the pdf dependence, leaving 
the pure BFKL effect. This turns out to be rather 
tricky [ |) , because the desired cancellations occur only 
at lowest order. Therefore we consider the ratio 

da(^/s 2 ,A 2 ) 

with A2 defined such that R12 — 1 in QCD lowest- 
order, the result is shown in Figure and we see that 
the kinematic constraints strongly affect the predicted 
behavior, not only quantitatively but sometimes qual- 
itatively as well. More details can be found in [ raj . 
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Figure 2. The ratio R\2 of the dijet cross sec- 
tions at the two collider energies y /s 1 = 630 GeV and 
y/s 2 = 1800 GeV, as defined in the text. The curves 
are: (i) the BFKL MC predictions (solid curve), (ii) 
the 'naive' BFKL prediction (dashed curve), and (iii) 
the asymptotic QCD leading-order prediction (dotted 
curve) R12 = 1. 
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3. CONCLUSIONS 

In summary, we have developed a BFKL Monte 
Carlo event generator that allows us to include the sub- 
leading effects such as kinematic constraints and run- 
ning of a s - We have applied this Monte Carlo to dijet 
production at large rapidity separation at the Teva- 
tron. We found that kinematic constraints, though 
nominally subleading, can be very important. In par- 
ticular they lead to suppression of gluon emission, 
which in turn suppresses some of the behavior that 
is considered to be characteristic of BFKL physics. It 
is clear therefore that reliable BFKL tests can only 
be performed using predictions that incorporate kine- 
matic constraints. 
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